Archive/NIMBioS.code/fastsimcoal/test hist ev check.R

In christianparobek/skeleSim: Genetic Simulation Engine

rm(list = ls())
source("fsc.set.params.R")


hist.ev <- rbind(
  c(2000, 1, 2, 0.05, 1, 0, 0),
  c(2980, 1, 1, 0, 0.04, 0, 0),
  c(3000, 1, 0, 1, 1, 0, 0),
  c(15000, 0, 2, 1, 3, 0, 0)
)

fsc.histEvConverges(hist.ev, c(20000, 5000, 10000), c(0, 0, 0))


hist.ev2 <- rbind(
  c(1000, 0, 0, 0, 1, 0, 1),
  c(10000, 1, 0, 1, 10, 0, 1)
)
fsc.histEvConverges(hist.ev2, c(1000, 1000), 0)


hist.ev3 <- rbind(
  c(70000, 0, 1, 1, 1, 0, 0)
)
fsc.histEvConverges(hist.ev3, c(2000, 1), 0)


hist.ev4 <- rbind(
  c(100000, 1, 0, 1, 1, 0, 1),
  c(80000, 2, 1, 1, 1, 0, 1),
  c(10000, 3, 2, 1, 1, 0, 1),
  c(50000, 3, 1, 1, 1, 0, 0)
)
fsc.histEvConverges(hist.ev4, rep(10000, 4), 0)


hist.ev5 <- rbind(
  c(100000, 1, 0, 1, 1, 0, 1),
  c(80000, 2, 1, 1, 1, 0, 1)
)
fsc.histEvConverges(hist.ev5, rep(10000, 4), 0)


hist.ev6 <- rbind(
  c(10000, 8, 0, 0.5, 1, 0, 0),
  c(5000, 0, 1, 0.05, 1, 0, 1),
  c(5000, 8, 1, 0.05, 1, 0, 1),
  c(5000, 0, 2, 0.05, 1, 0, 1),
  c(5000, 8, 2, 0.05, 1, 0, 1),
  c(5000, 0, 3, 0.05, 1, 0, 1),
  c(5000, 8, 3, 0.05, 1, 0, 1),
  c(5000, 0, 4, 0.05, 1, 0, 1),
  c(5000, 8, 4, 0.05, 1, 0, 1),
  c(5000, 0, 5, 0.05, 1, 0, 1),
  c(5000, 8, 5, 0.05, 1, 0, 1),
  c(5000, 0, 6, 0.05, 1, 0, 1),
  c(5000, 0, 7, 0.05, 1, 0, 1)
)
fsc.histEvConverges(hist.ev6, c(200000, rep(2000, 7), 50000), 0)